Infinitely Extended Kac Table 
of Solvable Critical Dense Polymers 



Or 



Paul A. Pearce 

Department of Mathematics and Statistics, University of Melbourne 
Parkville, Victoria 3010, Australia 

J0rgen Rasmussen 

School of Mathematics and Physics, University of Queensland 
St Lucia, Brisbane, Queensland 4 072, Australia 

Simon P. Villani 

Department of Mathematics and Statistics, University of Melbourne 
Parkville, Victoria 3010, Australia 

P.Pearce@ms.unimelb.edu.au, J.Rasmussen@uq.edu.au 
S.Villani@ms.unimelb.edu.au 



Abstract 

Solvable critical dense polymers is a Yang-Baxter integrable model of polymers on the square 
lattice. It is the first member CA4(1, 2) of the family of logarithmic minimal models CAi(p,p'). 
The associated logarithmic conformal field theory admits an infinite family of Kac represen- 
tations labelled by the Kac labels r,s = 1,2,.... In this paper, we explicitly construct the 
conjugate boundary conditions on the strip. The boundary operators are labelled by the Kac 
fusion labels (r, s) = (r, 1) ® (1, s) and involve a boundary field £. Tuning the field £ appropri- 
ately, we solve exactly for the transfer matrix eigenvalues on arbitrary finite-width strips and 
obtain the conformal spectra using the Euler-Maclaurin formula. The key to the solution is an 
inversion identity satisfied by the commuting double-row transfer matrices. The transfer matrix 
eigenvalues are classified by the physical combinatorics of the patterns of zeros in the complex 
spectral-parameter plane. This yields selection rules for the physically relevant solutions to the 
inversion identity which takes the form of a decomposition into irreducible blocks corresponding 
combinatorially to finitized characters given by generalized g-Catalan polynomials. This decom- 
position is in accord with the decomposition of the Kac characters into irreducible characters. 
In the scaling limit, we confirm the central charge c = —2 and the Kac formula for the conformal 
weights A r s = ( 2r ~^ ~ 1 for r, s — 1, 2, 3, . . . in the infinitely extended Kac table. 
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1 Introduction 



Materials such as plastics, nylon, polyester and plexiglass are made from polymers [1, 2, 3] 
which consist of very long chain molecules with a large number of repeating structural units 
called monomers. Polymers exist in low- or high-temperature phases which are characterized as 
either dense or dilute. Polymers are dense if they fill a finite (non-zero) fraction of the available 
volume in the thermodynamic limit. 

The modern era of two-dimensional polymer theory began in the late eighties [4, 5, 6] 
when Saleur and Duplantier initiated the study of polymers as a conformal field theory (CFT). 
From the viewpoint of lattice statistical mechanics, polymers are of interest as prototypical 
examples of systems involving (extended) non-local degrees of freedom. The non-local nature of 
these degrees of freedom has a profound effect on the associated conformal field theory (CFT) 
obtained in the continuum scaling limit. Indeed, the associated CFT is logarithmic [7] in the 
sense that certain representations of the dilatation Virasoro generator L Q are non-diagonalizable 
and exhibit Jordan cells. 

The first member £Ai(l,2) of the Yang-Baxter integrable family of logarithmic minimal 
models CAi(p, p') [8] on the square lattice is solvable critical dense polymers [9]. This model 
was solved exactly in [9] for certain boundary conditions of type (1, s) on a strip. The model has 
also been solved exactly on the cylinder [10] and its integrals of motion and Baxter Q-operator 
studied in [11]. With appropriate W-extended boundary conditions on the strip [12], the model 
also gives a lattice realization of symplectic fermions [13]. We refer to the bibliographies of these 
papers for related approaches to the problem and for more detailed references to the relevant 
literature. 

The conformal Virasoro algebra of solvable critical dense polymers is known [8, 9] to admit 
Kac representations labelled by Kac labels (r, s) in an infinitely extended Kac table 

Kac table: (r, s) : r, s = 1, 2, 3, . . . (1.1) 

as shown in Figure 1. From the fusion rule (r, s) = (r, 1) £g> (1, s), the finitized partition functions 
on a strip are 

43)Kr,)fo) = 4$i<i,i)fo) = *o$i<i,.)fo) = 4!ki)(?) = xSH?) (i-2) 

where Xr N s\q) are the finitized characters of the Kac representations. In the thermodynamic 
limit, the conformal partition functions are given by Virasoro characters 

%,i)l(rv0(?) = X^\q) = XrAl), r,s = 1,2,3,... (1.3) 

In this paper, we explicitly construct the conjugate boundary conditions on the strip labelled 
by (r, s) and solve for the conformal spectra given by X r ,s(q)- This is achieved by studying 
the model on a finite-width strip subject to the (r,s) = (r, 1) <g> (l,s) boundary conditions 
which involve a boundary field £. Remarkably if, following [14], the boundary field is tuned 
to the special value £ = A/2 = tt/4, the normalized double-row transfer matrices satisfy a 
functional equation in the form of a single inversion identity which is universal in the sense that 
it is independent of (r, s). This enables us to obtain the exact properties of the model on a 
finite lattice. The conformal properties are then readily accessible from finite-size corrections. 
In particular, in the continuum scaling limit, we confirm the central charge c = —2 and the 
conformal weights in the infinitely extended Kac table 

(2r — s) 2 — 1 

Ar,*=- ^ , r,s = 1,2,3,... (1.4) 
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The layout of this paper is as follows. In Section 2, we define the solvable model of critical 
dense polymers, the planar Temperley-Lieb (TL) algebra [15, 16] on which it is built and the 
double-row transfer matrices. We also define vector spaces of link states on which the transfer 
matrices act and relate these to the (r, s) boundary conditions. The boundary opperators are 
constructed in Appendix A. In Section 3, we review the c = — 2 CFT associated with critical 
dense polymers including a description of the Kac representations labelled by (r, s) and the 
decomposition of their characters into characters of irreducible representations. In Section 4, we 
present the inversion identity. Its proof in the planar TL algebra is relegated to Appendix B. 
We solve exactly the inversion identity on finite width strips to obtain the general solution for 
the eigenvalues of the transfer matrices. The physically relevant solutions, in the (r, s) sectors, 
are obtained empirically by selection rules encoded in irreducible blocks consisting of g-Catalan 
polynomials realized as suitable combinations of double-column diagrams. Combinatorics of the 
related g-Narayana numbers are given in Appendix C. The finitized (r, s) conformal characters 
are obtained explicitly from finite-width partition functions on the strip. The identities used to 
simplify the g-Catalan expression are given in Appendix D. It is also verified that the finitized 
Kac characters yield the quasi-rational characters of [8] in the continuum scaling limit. In 
Section 5, the finite-size corrections to the eigenvalues are obtained by applying the Euler- 
Maclaurin formula yielding the bulk and boundary free energies and confirming the central 
charge c = —2 and conformal weights A r s = ( 2r ~^ -1 , Lastly, we obtain the Hamiltonian limit 
of the double-row transfer matrices in the planar TL algebra and give its spectra. We conclude 
with a brief discussion in Section 6. 



2 Lattice Model 

2.1 Critical dense polymers 

Solvable critical dense polymers on the square lattice describes the statistics of extended non- 
local degrees of freedom or connectivities representing the polymers. Since the polymers are 
very long, the degrees of freedom are not allowed to form closed local loops. In the language of 
the logarithmic minimal models [8], the loop fugacity vanishes /3 = 2 cos A = with the crossing 
parameter A = n/2. In accord with being in the dense phase, the polymers or non-local degrees 
of freedom visit every face of the lattice exactly twice. More explicitly, consider a finite strip 
consisting of N columns and M double rows of faces. An elementary face of the lattice can 
assume one of two configurations with different statistical weights 




(2.1) 



where the arcs represent local segments of polymers. A typical configuration looks like 
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Once the boundary conditions on the left and right edges have been specified, one is left with 
non-local degrees of freedom corresponding to a number of polymers on a strip. 



2.2 Planar Temper ley- Lieb algebra 

Algebraically, solvable critical dense polymers is described by the planar Temperley-Lieb (TL) 
algebra [15, 16]. An elementary face of the square lattice is assigned a face weight according to 
the configuration of the face. The two possible configurations with their associated weights are 
combined into a single face operator as 



u 



sin w 



+ COSM 



(2.3) 



where the lower left corner has been marked to fix the orientation of the square. 

The planar TL algebra is a diagrammatic algebra built up from elementary faces and 
triangles. The triangles enter through boundary conditions on the strip. The faces and triangles 
are connected such that the midpoint of an outer edge of a face, called a node, can be linked 
to a node of any other (or even the same) face as long as the total set of links make up a non- 
intersecting planar web of connections. Two relevant local properties of the planar TL algebra 
are the inversion relation 




cos 2 f 




(2.4) 



and the Yang-Baxter equation (YBE) [17] 
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(2.5) 



These are identities for 2- and 3-tangles, respectively, where a fc-tangle is an arrangement of 
faces and triangles with 2k free nodes. 



2.3 Double-row transfer matrix 

In the planar algebra, the double- row transfer matrix D(u) on a strip with (r, s) boundary 
conditions is defined as the (N+p+s — 2)-tangle with N bulk columns and p + s — 2 boundary 
columns where p, s > 1 



D(u) 
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(2.6) 



N 



p-1 



Two boundary seams of type (r, 1) and type (l,s) on the right are delimited from the bulk 
on the left and from each other by dashed vertical lines. By construction (Appendix A), the 
boundary operator 



(r,s) 



\ \ V'"\ \ s 



(M) 



®(M) 



(2.7) 



p-1 

satisfies the Boundary Yang-Baxter Equation (BYBE) 



s - 1 





(2- 



The bulk consists of N columns with parity fixed by iV = p + s mod 2. The r-type seam consists 
of p — 1 columns and the s-type seam consists of s — 1 columns. All faces in the bulk and r-type 
seam have the standard orientation with the bottom-left corner marked and spectral parameters 
as indicated. The column inhomogeneities are 



6 = e+jA 



(2.9) 



where £ is a boundary field. The face configurations in the s-type seam are fixed as shown. 
An s-type seam is just an r-type seam in the limit £ — > ioo. For critical dense polymers, the 
crossing parameter is A = |. In this case, the non-negative integer p can be chosen with either 
even or odd parity and is related to r by 



P 



2r, p even 
2r — 1, p odd 



(2.10) 



For the special value of the boundary field £ = j, in particular, this is the appropriate choice 
of p to realize the r-type conformal boundary conditions. More generally, for real £, the lattice 
boundary conditions can converge in the continuum scaling limit to conformal fixed points of 
different r-type boundary conditions depending on intervals of £ corresponding to respective 
basins of attraction. 

The double-row transfer matrix D(u) coincides with the double- row transfer matrices 
introduced in [8] for general fugacity = 2 cos A. In the planar algebra, multiplication of D(u) 
with D(v) means vertical concatenation of the two tangles. Using diagrammatic arguments 
based on [18], it follows that D(u) is crossing symmetric and gives rise to a commuting family 
of transfer matrices 



D(X-u) = D(u), [D(u),D(v)] = 



(2.11) 
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2.4 Link states and matrix representations 



Since the double-row transfer matrix is constructed as a planar tangle, a matrix representation 
is needed to obtain spectra. A matrix representation of the double- row transfer tangle D(u) is 
obtained by acting with the planar tangle D(u) from below on a suitable vector space of link 
states. We use the same notation for the tangle and the matrix representation. For the (r, s) 
boundary conditions, the link states must be restricted so that there are no half-arcs closing in 
the p — 1 columns of the r-type seam and no half-arcs closing in the s — 1 columns of the s-type 
seam. Half-arcs are however allowed to close between the r- and s-type seams. Let Vf£P denote 
the vector space of these link states with N + p + s — 2 nodes. Its dimension is given by 



N 

N-p+i 
2 



N 

N-p- 



N 

N-2r+s 



N 

N-2r+s+l 



N 

N-2r-s I ■ 



N 

N-2r-s+l 



P 



2r 



2r 



(2.12) 



(TV) 



Using dashed vertical lines to separate the bulk and boundary seams, the six link states for V 3 3 
with N = 4 are 




i « 'i 

i i 

i i 

i i 

i i 




i i 




(2.13) 



A link state on In nodes consists of n half-arcs. We call a half-arc connecting neighbouring 
nodes a small half-arc. If there is a half-arc anywhere in a link state, then either the half-arc 
is itself a small half-arc or there must be a small half-arc somewhere inside the given half-arc. 
Within r- and s-type seams, small half-arcs have a push-through property. Specifically, if there 
is a half-arc (and therefore a small half-arc) closing inside the top of an r-type seam then, by 
a simple trigonometric identity, there must also be a small half-arc in the out link state at the 
bottom of the seam 



v+X 



s(A - v)s(X + v) 




+ s(\-v)s(- 




(2.14) 



+ 0s(v)s(-v) 



+ s(v)s(v + A) 




s(A - v)s(X + v) 




This holds for arbitrary v,X with (3 = 2 cos A and s(v) = sin v. The same is trivially true 
for s-type seams. We want the double-row transfer matrices to act from the space Vj£P back 
to itself. It follows from the push-through property that, if we restrict the out states at the 



bottom to Vp\P , then the in-states at the top must also be restricted to Vp" s ' and likewise for all 
intermediate link states. This restriction preventing half-arcs closing within a seam is equivalent 
to the action on the seam of the Wenzl- Jones projectors [19, 20] when they exist. But the half- 
arc restriction makes sense even when the Wenzl- Jones projectors do not exist which is the case 
for critical dense polymers. 



(TV) 
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3 Conformal Field Theory 
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Figure 1: Kac table of conformal weights of critical dense polymers. Each position (r, s) of the 
table is associated with a Kac representation, a conformal dimension A r8 and an associated 
Kac character Xr,s(o)- The (r, s) representations indicated with a red quadrant are irreducible 
representations with characters X r ,s(q) — ch r s (g). The irreducible representations at positions 
(1, 2k) = (Jfe, 2), k = 2, 3, . . . are identified. 



3.1 Spectrum and characters of critical dense polymers 

In the continuum scaling limit, the finite-size spectra of the transfer matrices of critical dense 
polymers are described [8, 9] by a logarithmic CFT with central charge c = —2 and conformal 
weights A r s given by (1.4). This CFT is non-unitary and non-rational. 

In the Virasoro picture, the CFT data is neatly encoded in the infinitely extended Kac table 
shown in Figure 1. The (r, s) conformal characters 

s_ — +A 00 

xrM) = fe^v-* - qAr ^ = q ~i^- 1 ^ - ^ = ri( i - qk ) 

\H)oo \H)oo k=l 

are the generating functions for the conformal finite-size spectra of the double-row transfer 
matrices on the strip with (r, s) boundary conditions. The associated representations are called 
Kac representations and described in Section 3.2. We reserve the notation ch rs (q) = X r , s (q) 
with s — 1, 2 to denote the characters of the irreducible representations in the first two rows of 
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the Kac table 

ch 1 ,.(q) = ^ l-g rs , s = l,2 (3.2) 

Woo 

The set of Kac representations is not exhaustive. Although we do not do it here, critical 
dense polymers can also be considered in the W-extended picture with a W-extended conformal 
symmetry [21, 22, 23]. In this picture, it suffices to consider a finite number of representations 
and extended boundary conditions which respect the enlarged symmetry algebra. In fact, 
critical dense polymers with W-extended boundary conditions [24] gives a lattice realization of 
symplectic fermions [13]. The structure of the corresponding W-extended Kac representations 
is discussed in [25]. 



3.2 Description of Kac representations 

The Virasoro character of the Kac representation (r, s) is given by (3.1). This is recognized as 
the character of the quotient module 

Qr,s = V r ,./Vr,-„ T, S G N (3.3) 

obtained from the highest-weight Verma module V r>a of conformal weight A rjS by quotienting 
out the submodule V r - S at level rs. However, this does not imply that the Kac representation 
(r, s) can be identified with the quotient module Q r , s - Indeed, a characterization of the Kac 
representations as modules over the Virasoro algebra was recently proposed in [26], where 
it was argued that they can be understood as finitely generated submodules of Feigin-Fuchs 
modules [27]. A necessary condition for this conjecture to be true is that the characters of these 
finitely generated submodules match the characters (3.1) of the corresponding quotient modules. 
As seen below, this is indeed the case. 

For s even, the Kac representation (r, s) = (r, 2k) is fully reducible [28] 

r+Jfe— 1 

(r,2k) = M h2 (3.4) 

j=\r-k +l,by 2 

where M m>n is the irreducible highest-weight module of conformal weight A mn , and the 
corresponding Virasoro character is given by 

r+fc-l 

X r ,2k{q) = ^2 cn i.2(<?) (3-5) 

J=|r-fc|+l,by2 

For s odd, on the other hand, the Kac representation (r, s) = (r, 2k + 1) is believed [26] to be a 
reducible but indecomposable rank 1 representation characterized by 

(r, s) : ->■ M fc _ r+2 ,i <- M fc _ r+3i i M k+r -i,i -> M k + r ,i, 2r < s 

(3.6) 

(r, s) : M r _ M <- M r - k+1A -)■ M r _ fc+2 ,i M r+k -i,i ->■ M r+kjl , 2r > s 

The corresponding Virasoro character is given by 

r+k 

x r ,2k+i(q) = ch jM ( 3 - 7 ) 



s 



Using the explicit expressions (3.2) for the irreducible characters ch rtS (q), and observing the 
telescoping property of their sums in (3.5) and (3.7), it is readily verified that these sums are 
equivalent to the general expression (3.1) for the Kac characters X r ,s(o)- For our purposes here, 
we write the decomposition of Kac characters into irreducible characters as 



( i 
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XrM 



(2r+s-l) Min[2r 



c1i m(<?) = Ch |(|2r-s|-l+2fc),l(?) = Cll |(2r-s-l+2fc),l(^)' S ° M 

fc=±(|2r-s|+l) k=l k=l 

(3.1 

i(2r+s-2) Min[r,s/2] s/2 

Y ch fcj2 (g) = ^ chi ( | 2r _ s |_ 2+4fe)2 (g) = ^chi (2r _ s _ 2+4fc)2 (g), seven 

fc=i(|2r-s|+2) k=l k=l 



The last form of these identities uses the extended definitions 

ch 0iS (g) = 0, ch r>s (q) = -ch_ r>a (g), r < 0, s = 1,2 (3.9) 

It is noted that the set of irreducible Kac representations is given by 

(r,l), (r,2), (1,2*;), r,kN (3.10) 

However, since an irreducible representation is uniquely characterized by its associated conformal 
weight, we identify the Kac representations 

(1,2k) = (k,2), keN (3.11) 

The set of inequivalent irreducible Kac representations is thus associated with the first two rows 
of the extended Kac table, that is, s = 1, 2 in (r, s). 

4 Exact Solution on a Finite Strip 

4.1 Inversion identity 

Remarkably, the double-row transfer matrix satisfies an inversion identity in the planar TL 
algebra. This inversion identity is unique to critical dense polymers among the logarithmic 
minimal models. 

Inversion Identity: For A = | ; the planar tangle D(u) (2.6) satisfies the inversion identity 
D(u)D(u + A) = - sin 2 2u ^\u, £) V ^(u + A, £) f p (u, f) I (4.1) 

where the right side is an entire function of u, I is the vertical identity tangle, 

p-i 

V W (u, = J] sin( M + sin( M - (4.2) 



and 



,'cos 2Ar u — sin 2JV M sin(p — l)Acos 2Jv w 
/p(«>0 = 5 + 



cos 2 u — sin u cos(u + £) cos(m — £ p+1 ) 

/cos 2Ar w — sin 27V -u ^ sin(p — l)Asin 2Ar -u . ^ 
cos 2 u — sin 2 u sin(w + £) sin(u — £ p +i) 
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This inversion identity is proved in the planar algebra in Appendix B. For p odd, the 
function f p (u,£) reduces to 



2N ■ 2N \ 2 

cos u — sin u \ 
cos 2 u — sin 2 u / 



p = 2r - 1 



while for p even, we have 

[(cos 2N u — sin 2N u) cos 2£ — (cos 2N u + sin 2N u) cos 2u\ ' 



f P (u,0 



cos 2 2u (cos 2 2£ — cos 2 2u) 
For the special value of £, 

^( W> f) = (icos2«) 
and the general inversion identity (4.1) reduces to 



p = 2r 



p-l _ A _ 7T 



2 4 



cos^ M-sm M 



4.2 Exact solution for the eigenvalues 

Let us introduce normalized transfer matrices by 

' 2P- 1 D(u) 
sin 2u cos p ~ 2 2u ' 

d{u) 

2' J - 1 D(u) 



p = 2r 



p = 2r-\ 



y sin 2u cos p 1 2u ' 

The eigenvalues d(u) of the normalized double-row transfer matrices d(u) then satisfy 

(cos 2iV M + sin 2Af M ) 2 , 



d{u)d{u + A) 



2N • 2N \ 2 

cos u — sin u 
cos 2 u — sin 2 w 



p = 2r 
p = 2r - 1 



subject to the initial condition and crossing symmetry 

d(0) = 1, d(A - it) = 



(4-4) 



(4.5) 



(4.6) 



D{u)D{u + A) = - sin 2 2 M (\ cos 2 M ) 2p " 2 f ^gg^-J- ^^^ ^1, f) G N (1.7) 



(4.t 



(4.9) 



(4.10) 



We observe that the form of this inversion identity for p odd exactly agrees with the inversion 
identity solved in [9] for r = 1. The p-dependent order-1 normalization terms in (4.8) exactly 
cancel the additional boundary free energy arising from the r-type seam. 
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The finite-size analysis of (4.9) depends on the following factorizations 
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cos u — sin w 
cos 2 w — sin 2 w 



--1 



]'[ (1 -sin 2 f sin 2 2u), 



p = 2r, s odd 
p = 2r, s even 

p — 2r — 1, s odd 



(4.11) 



jV-l 
2 

n 



(4.12) 



sm 



2 (2j-l)7, 



2N 



sm' 



2m), p = 2r-l, 



s even 



It is now straightforward to solve the functional relation (4.9), subject to (4.10), by sharing out 
the zeros. Recalling that iV + p + s is even, we find 



d{u) 



jQ (l + £j sin^ sin 2u) (l + pj sin^ sin 2uj . 



s odd 



N_ 

2 - 



LfJ 



(4.13) 



n(!+ 



• (2jf — 1)tt 
e 3 SmLJ 2N 



sin 2m) (l + /ij sin ^ J 2jV sin 2m) , s even 



where e 2 = u 2 = 1 for all j. The appearance of two sets of parameters {ej} and {u^} stems from 
the overall squaring in (4.9) and allows for double zeros. 

For either parity of s in (4.13), the maximum eigenvalue is obtained for 6j = pj = 1 for all 
j and corresponds to the groundstate in the (r, s) = (1,1) or (1,2) sector. Excited states are 
generated by switching a (finite) number of the parameters ej, pj from +1 to —1. The number of 
possible excitations following from (4.13) is 2 2 LC^— x ) / 2 J anc l 2L iV / 2 l, respectively, compared with 
the number of link states in a given p, s sector given by (2.12). The inversion identity fixes the 
form of the eigenvalues but not their multiplicities. We thus need a set of selection rules to 
determine the groundstate and excitations, along with their multiplicities, that actually appear 
in the spectrum in a given sector. This information is neatly encoded by specifying the patterns 
of complex zeros {ej,pj} for the allowed eigenvalues d{u). 

From the crossing symmetry and periodicity, the zeros Uj of an eigenvalue d(u) come in 
complex conjugate pairs in the complex M-plane and appear with a periodicity tt in the real part 
of u. From (4.13) 



TT 



, - (2 + Uj)- ± -In tan— + ttZ 
4 2 2 



(4.14) 



where Uj is €j or pj, and 



22L 
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(2j-l)7T 

2N ' 



s odd 



s even 



(4.15) 
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A typical pattern for N = 12 is 
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(4.16) 



where 



y.i 



% tj 

■— m tan — 
2 2 



(4.17) 



The analyticity strip — j < Keu < ^ is called the physical strip. All zeros in the strip lie either 
on the boundary or on the vertical centre line (1-strings). In contrast to the usual situation in 
the context of Bethe ansatz, these zeros can occur either as single or double zeros. A single zero 
is indicated by a grey dot while a double zero is indicated by a black dot. Counting a double 
zero twice, the number of zeros with fixed imaginary value ±yj and real part either | or ^ is 
two. This follows straightforwardly from (4.14) since Uj can be either ej or fij, as depicted in 
(4.16). 

From complex conjugation symmetry, it follows that the full pattern of zeros is encoded in 
the distribution of 1-strings in the lower half-plane. This distribution is actually a sum of two 
- one governed by excitations administered by e and one by fi. The separation of the 1-strings 
into contributions coming from e and /i, respectively, is superfluous but nonetheless a helpful 
refinement later for the description of the selection rules. It is illustrated here 



i = 2 
i = i 



< — y 



(4.18) 



where the left side corresponds to (4.16) while the two columns to the right encode separately 
the e and /x excitations. Not all such double-column configurations will appear in a given sector. 
A characterization of a set of admissible two-column configurations provides a description of 
the selection rules. This combinatorial designation of the physical states is termed 'physical 
combinatories'. As we will see, the double-column configurations provide a natural basis for 
defining the associated finitized characters. These double-column configurations and finitized 
characters also arise in the context of the dimer model [29]. 
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4.3 Double-columns and (/-Catalan polynomials 



A single- column configuration of height M consists of M sites arranged as a column. The 
sites are labelled from the bottom by the integers 1,2,..., M. Each site can be occupied or 
unoccupied. The configuration is specified by listing the labels of the occupied sites in descending 
order. A double-column configuration S of height M consists of a pair (L, R) of single-column 
configurations of height M. For the double-column configuration on the right side of (4.18) 



S 



(L,R), 



(3), 



R 



(4,3,1) 



(4.19) 



If m = \L\ = §$^=1(1 ~~ e i) an d n = \R\ = |X)i=i(l ~ A*j) indicate the number of occupied 
sites in the left and right columns, respectively, the energy E = E(S) of a double-column 
configuration is given by 



E(S) = E(L,R) 



where we will ultimately take 



Lj — Rj — Ej 



1 

2' 



j is occupied with s odd 
j is occupied with s even 
j is unoccupied 



In Unitized characters, we will also associate the monomial 

q E(S) 



(4.20) 



(4.21) 



(4.22) 



to the double-column configuration S. The formal parameter q will ultimately be the modular 
nome. 

A double-column configuration S = (L, R) is called admissible if L ^ R with respect to the 
partial ordering 

L r< R if Lj < Rj, j = l,...,m (4.23) 



which presupposes that 



< m < n < M 

Admissibility can be simply characterized geometrically as illustrated in the example 



(4.24) 




(4.25) 



One draws line segments between the occupied sites of greatest height in the two columns, 
then between the occupied sites of second-to-greatest height and so on. The double-column 
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configuration is now admissible if it does not involve line segments with a strictly negative slope. 
Thus, in an admissible double-column configuration, there are either no line segments (m = 0) 
or each line segment appears with a non-negative slope. 

Combinatorially, we define the (generalized) g-Narayana numbers ^ ^ ^ as the sum of the 

monomials associated to all admissible double-column configurations of height M with exactly 
m and n occupied sites in the left and right columns respectively 



M 
m,n/ q 



E 

S:\L\=m,\R\=n 



E(S) 



(4.26) 



In Appendix C, we show that these polynomials admit the closed-form expressions 

,±m(m+l)+|n(n+l) 



M 

m,n/ g 



q2' 





~ M~ 




~M~ 




M 




M 




( 


m 


q 


n 


_ qn—m+l 
1 


m — 1 


9 


n + 1 


.) 



(4.27) 



where 



M 
rn 



is a g-binomial (Gaussian polynomial). These coincide with g-Narayana num- 
bers [30, 31] when m = n. 

Following [10], for r > 0, we define the following (generalized) g-analog of Catalan numbers 



M-r+l 



c 



mA 



(?) 

C' Mt M = q 



E 

m=0 



M 

m, m+r- 



r(r-l) (1 — q T ) 



;i - q M+1 ) 



2M + 2 
M + 1 - r 



s odd 



M-r+l 

E «■ 



M 



q 



(r-D 2 (1 - q 2r ) 



(4.28) 



m=0 



2M+ 1 
M+ 1 - r 



6' even 



,m,m+r-l/g " (1 - g^+r+i) 

It is convenient to extend these definitions using symmetry 

^o(g) = C' M>0 (q) = 0, C M , r (g) = -C M ^ r (q), C' M;r (q) = -C' M> _ r (q), r<0 (4.29) 

The polynomials CM,i{q) coincide with one common definition of g-Catalan numbers. In the 
sequel, we will just refer to the polynomials in (4.28) as g-Catalan polynomials. Although we 
refer to these as g-Catalan polynomials we stress that, as defined for r even, C' Mr (q) is actually 
a polynomial in y/q due to the prefactor. The g-Catalan polynomials (4.28) are the spectrum 
generating functions of the irreducible building blocks for the finitized characters in the case 
that s is odd or even respectively. Despite the minus signs in (4.27) and (4.28), the g-Narayana 
numbers and g-Catalan polynomials are fermionic in the sense that they only have non-negative 
coefficients. 

Combinatorially, for r > 1, the (generalized) g-Catalan polynomial CM,r{q) is the sum 
of the monomials associated to all admissible double-column configurations of height M with 
n — m = r — 1, that is, exactly r — 1 more occupied sites in the right column compared to the 
left column. With the energy assignment Ej = j appropriate for s odd, we have 



CM,r{q) 



E " 

S: \R\-\L\=r-l 



E(S) 



(4.30) 



Similarly, for r > 1, the (generalized) g-Catalan polynomial C' Mr (q) has precisely the same 
combinatorial interpretation but with the energy assignment Ej = 3 — \ appropriate for s even. 
The case M = r = 2 is illustrated in Figure 2. 
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Figure 2: Illustration of the combinatorial interpretation of the g-Catalan polynomials C^ifi) = 
q + q 2 + q 4 + q 5 and C 2 2 (q) = + q + q 2 + q 3 ) ior M = r = 2. The energy assignment 

is Ej = j,j — | as appropriate for s odd, s even respectively. The sum is over all admissible 
double-column configurations of height M with r — 1 more occupied sites in the right column. 



The g-Catalan polynomials are simply related to finitized irreducible characters 

= q~^C M , r (q)i ch^ 2 I] (q) = q~&~sC' M>r (q) 

Using the result 



lim 

M->oo 



M 

m 



(q\ 



(q) m = l[(i-q k 



it follows from (4.28) that in the thermodynamic limit 

lim ch^ f) (g) = ch nl (g), lim ch^\q) = c\ 2 (q) 



(4.31) 



(4.32) 



(4.33) 



4.4 Selection rules 

We now turn to the description of the selection rules. They depend on the parities of p and 
s with N + p + s even. The relationship between p and r is as given in (2.10). Based on a 
body of empirical data obtained by examining patterns of zeros by numerically diagonalizing 
the double- row transfer matrices out to sizes N = 14 with 1 < p, s < 8, we make the following 
conjecture. 

Selection rules: For a system with N bulk columns, p — 1 r-type and s — 1 s-type boundary 
columns, the selection rules in the (r, s) sectors are equivalent to singling out the following sets 
of admissible double- column configurations in terms of the irreducible building blocks given by 
q- Catalan polynomials 



sodd: Z^J )l[r;S) (q) 



A 1 p-s-l+2k (q) , 

2 



q 24 22 

k=l 
s 

g _ 24 ^ [CV-2 p-s+2k (q) + q^C N -2 p-s-2+2 



fc=l 



seven: Z^^iq) 



c__ 1 

" 24 8 



_£__ 1 

" 24 8 



s/2 

k=l 

s/2 



JV p-s-2+4fc \L{ )i 
2 ' : 



p = 2r 
p = 2r-l 

p = 2r 



(4.34) 



(4.35) 



[C 'n-i P s-i+4k (q) + q 1 2C ' N ^ 1 e - £ -3+4fc (<?)] , p = 2r-l 



k=l 
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If s > 2r, the g-Catalan polynomials CM,r{q) and C' Mr (q) with r < explicitly cancel out in 
these sums against the same polynomial with r > 0. Although we do not do so, the sums in 
these cases can be restricted to the values of k that do not cancel out. In the first sum with 
s > 2r odd and p even, for example, the lower terminal could be replaced with k — s — 2r + 1. 

Notice that the excitations form conformal towers above the groundstate with integer 
spacing. To see this, let us define the excitation energies 



s odd 
s even 



and the modular nome q 



q 



-2-ET 



T 



M 

— sin 2u 

N 



(4.36) 



(4.37) 



where M/N is the aspect ratio and $ = 2u is the anisotropy angle [32] related to the geometry 
of the lattice. Then, for s odd, taking the ratio of (4.13) with precisely one €j or pj = — 1 to 
(4.13) with all €j = pj = +1 and taking the limit M, N — > oo with a fixed aspect ratio 5 = M/N 
gives 



lim 



1 



sin ^ sin 2u 



m,n-+oc \ l + sin ] -§ sin 2u 



M 



exp[— 2j7r 5 sin 2-u] = q 1 



J 



(4.38) 



Similarly, for s even, taking the ratio of (4.13) with precisely one 6j or fij = — 1 to (4.13) with 
fij = +1 and taking the limit M, N — > oo with a fixed aspect ratio 5 = M/N gives 



all €j 



1 



sin 



lim 

m,tv— »oo Vl + sin 



(2j-l)w N M 
2N 



(2j-l)7T 

2N 



exp[— (2j — 1)tc 5 sin2-u] = q 1 



E; 



J 



(4.39) 



The groundstates in the various sectors are given by the zero patterns associated with the 
minimum surviving energy configurations in (4.34) and (4.35) as shown in Figure 3. For 2r > s, 
this occurs in the g-Catalan polynomial with k — 1. In all cases, for p even or odd, it occurs in 
the leading g-Catalan polynomials CM,n-m+i(q) or C' M n _ m+1 (q) where the excess occupation in 
the right column over the left column is given by 



n — m 



\R\ - \L\ 



§(|2r-s| -1), sodd 
\{\2r — s\), s even 

The conformal weights associated to the minimal configurations are 



(4.40) 



An-m+1,1) S odd 

A n _ m+ i j2 , s even 



(4.41) 



Applying the identities proved in Appendix D, the g-Catalan decompositions (4.34) and 
(4.35) reduce to simple finitized characters 



{i,i)\(r, s )(i) ~ xiy{q) 



N 

N-2r+s 



N 

N-2r+s+l 



N 

N-2r~s 



N 

N-2r-s+l 



p = 2r 
p = 2r-\ 



(4.42) 
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Figure 3: Groundstate double-column configurations arranged in a Kac table for r, s = 1, 2, 3, 4. 
The continuation of the pattern for larger values of r and s is clear. The solid grey dots represent 
single 1-strings in the center of the analyticity strip. The excess occupation in the right column 
over the left column is given by (4.40) and the conformal weights by (4.41). There are no double 
zeros in the center of the analyticity strip for these groundstates. 



The results with p = r = 1 and s arbitrary exactly agree with the results of [9]. We also observe 
that the limit q — » 1 gives the correct counting of states (2.12). Taking the thermodynamic 
limit and again using (4.32), we obtain the desired result 

Z { i, mr ,s)(<l) = l r im =Xr,,(g), r,s = 1,2,3,... (4.43) 

Similarly, taking the thermodynamic limit directly on (4.34) and (4.35) using (4.31) and (4.33) 
gives the irreducible decomposition (3.8) of the Kac character Xr,s{o)- While these results 
depend on the conjectured selection rules, we point out that these selection rules were previously 
conjectured [9] for the case p — r — 1 with s arbitrary and subsequently proved by Morin- 
Duchesne [33]. 

5 Finite-Size Corrections 
5.1 Euler-Maclaurin 

The partition function of critical dense polymers on a lattice of iV columns and M double rows 
is defined by 

Z N , M = Trd(u) M = Y, d n{u) M = J2e~ ME " (u) (5.1) 
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where the sum is over all eigenvalues of the normalized transfer matrix d(u) including possible 
multiplicities and E n {u) is the energy associated to the eigenvalue d n {u). Although we suppress 
the dependence on (r, s), the form of the partition function (5.1) holds in all sectors. Conformal 
invariance of the model in the continuum scaling limit dictates [34, 35] that the leading finite-size 
corrections for large iV in the (r, s) sector are of the form 

2tx sin 2u / c \ 
E n (u) = -\nd n (u) ~ 2Nf bulk + f bdy + — (-— + A r>s + kj (5.2) 

Here f bu ik is the bulk free energy per face and f b d y is the boundary free energy (excluding 
the additional boundary free energy of the r-type seam which was explicitly removed by 
normalization of d{u)). The central charge of the CFT is c while the spectrum of conformal 
weights is given by the possible values of A, r s with excitations or descendants labelled by the 
non-negative integers k. As the analysis below will confirm, the asymptotic behaviour of the 
eigenvalues (4.13) are in accordance with (5.2). 
Let us introduce the function 

F(t) = ln(l + sin t sin 2u) (5.3) 

where for simplicity the u dependence has been suppressed. As a partial evaluation of the 
energies (5.2), we find 

In JJ (1 + ej sin tj sin 2u) ~ F(tj) - 2 sin 2u ^ tj (5.4) 
3 3 je£n 

where S n is the subset of j indices for which ej — — 1 and E n remains finite as N — > oo. Likewise, 
A4 n denotes the subset of j indices for which fij = —1 and it also remains finite as iV — > oo. 
The argument tj is defined in (4.15). The endpoint and midpoint Euler-Maclaurin formulas are 

m 1 f b 1 h 

J2 F (a + kh) = - F(t)dt+-[F(b) + F(a)} + -[F'(b)-F'(a)} + 0{h 2 ) 

k=o m Ja (5 5) 

Y J F{a+(k-\)h)= l 1 -f b F(t) dt - ^\F'{b) - F\a)} + 0(/> 2 ) 
k=i Ja 

where h = (b — a)/m. These formulas are valid since F(t) and its first two derivatives are 
continuous on the closed interval [a, b}. Applying these formulas, we find 



E n (u) ~ 2Nf bulk + f bdy 

V24 



2tt sin 2u 



' f— + j )' sodd 



j&£n j&Mn 



je£n jeM„ 
In these expressions, the bulk free energy per face is given by 



(5.6) 



s even 



i r 1 i r 1 / 1 \ 

fbuik = / ln(l + sint sin2u)dt = |ln2 / In h sin 2t* J dt (5.7) 

W 2 71 J Vsmt / 

whereas the boundary free energy is given by 



fi 



bdy 



(-2f Mk -F(z), p = 2r, sodd 

0, p = 2r, s even 

-4fbui k -F(f), p = 2r-l, sodd 

-2 f talk, p = 2r-l, seven 



(5i 
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with F(f) = ln(l + sin 2m). 

The link between the lattice model and the characters of the CFT is governed by the 
modular nome q defined by (4.37). 



5.2 Hamiltonian limit 



Since the normalized transfer tangle d(u) reduces to the identity tangle for u = 0, we define the 
Hamiltonian limit of d{u) as the tangle appearing as the first-order term in the Taylor series of 
d{u) centred at u = 0. We thus consider 



d{u) 



2uH + 0{u 2 



H 



1 d 

2 du 



In d(u) 



u=0 



(5.9) 



where the factor of 2 reflects that d(u) is a double-row transfer tangle. Employing techniques 
used in [8, 9], we find the Hamiltonian 



// 



N-l 
3=1 



3 (P) 
N 



(5.10) 



here written in terms of the linear TL algebra, with loop fugacity = 0, acting vertically on 
the vector space vj£P of link states with N + p + s — 2 nodes. The operator 



P 



(p) 



N 



( (p-2)/2 

E< 

fc=0 
(p-l)/2 

E( 



>(p-2fe-2)/2 p (2fc) 



p = 2r 



(5.11) 



(p_ 2 fc+l)/2 _(2fc-l) 
' e 7V ? 



2r-l 



fc=i 



where 



,(*) 



■ ■ ■ &N+k 



X\e N+J (5.12) 

is the generalized Temperley-Lieb projector discussed in Appendix A. The set of normalized 
eigenvectors of H is the same as that of d(u). While the eigenvalues of H and d(u) are 
different, the finite-size corrections are still given by (5.2) but with the bulk and boundary free 
energies replaced with their derivatives evaluated at u = 0, as in (5.9), and without the factor 
sin 2u in the conformal term. The conformal spectra are thus precisely the same. 



6 Conclusion 

In this paper, we have solved exactly a model of critical dense polymers on strips of arbitrary 
finite width. The spectra of this model have been obtained within the framework of Yang-Baxter 
integrabilty by solving exactly a functional equation in the form of an inversion identity satisfied 
by the double-row transfer matrices. The calculations have been carried out with the vacuum 
boundary condition (1, 1) on the left edge of the strip and for an infinite family of integrable 
boundary conditions, labelled by the Kac labels (r,s), on the right edge. Our study of the 
physical combinatorics, which is the classification of the physical states using the combinatorics 
of patterns of zeros in the complex w-plane, has revealed some interesting connections to 
g-Narayana numbers and g-Catalan polynomials. The conformal spectra have been obtained 
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from finite-size corrections using the Euler-Maclaurin formula. The results confirm the central 



charge c = —2, the conformal weights A f 
functions and Kac characters 



(2r-s) 2 -l 



and the form of the conformal partition 



q "H +Ar < s 



(?)c 



l-q r 



r,s = 1,2,3, 



(6.1) 



Importantly, this analysis fully establishes the infinitely extended Kac table of Figure 1 for 
CM (1,2) [8]. In closing we observe that, for (r,s) satisfying (2r — s) 2 < 8, the conformal 
dimensions also agree with the critical exponents (3 rjS = (2 — a)A r>s = A r s associated with 
generalized order parameters of the exactly solvable ^1,3 off-critical perturbation of dense 
polymers 
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A Construction and Expansion of Boundary Operators 

In this Appendix, we construct the r-type boundary operators K^ p '(u,^) for p — 1 columns. The 
construction is general giving the r-type boundary operators for all of the logarithmic minimal 
models CM.(p,p'). So, in particular, in this appendix the crossing parameter A is not restricted 
to 7r/2 and we use the notation 



s[u) 



S11111 

sin A ' 



p -p 

A = 7T 

pi 



(A.r 



The boundary operator consists of a planar diagram along with a restriction on the link 
states. The planar boundary diagram is 



K^(u,C) 



1 








-U-nfo 


v ip) (u,o _ 








u-ii 



(A.2) 



where there are p — 1 columns, the column inhomogeneities are £j = £ + jX and ^ is a boundary 
field. By construction, following [14] and using the push-through property (2.14), the boundary 
operator K^ p \u,^) is a solution to the boundary Yang-Baxter equation (2.8) for all p, u and £. 
A convenient normalization factor is 



p-1 



(A.3) 



As we will see, this ensures that K^(0,^) = I. 
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Opening up the boundary diagram gives 



7)(P)(u,£) 




(A.4) 



The boundary operator can therefore be written in the linear Temperley-Lieb algebra, acting 
from the bottom to the top, as 



r j ^\u^)K { p \u^) = X (u-^ 1 )X l {u-^ 2 )- ■ ■X^(u-^)X f ^(uHi)X ( ^(u+^) . ..X (u+^i) 

(A.5) 

where the face operators are 



Xj{u) = s(\ — u)I + s{u)ej 



(A.6) 



The index in Kq P \u } ^) indicates that the leftmost strand is in position j = 0. Note that a 
closed half- arc is allowed between the nodes in position and 1. 

The link states are restricted such that closed half-arcs between the rightmost p — 1 nodes 
are disallowed. Since this restriction applies to all link states it implies that, in expanding out 
the planar diagram, we can ignore connectivity diagrams containing half-arcs between the p — 1 
nodes on the lower edge of (A. 2). We thus have the diagram expansion 




The first diagram acts as the identity on the p — 1 vertical strands. Opening up the other 
diagrams, we see that they coincide with the TL words e ei . . . for k = 0, 1, . . . , p — 2. All 
other TL words are killed by the restriction on link states. The coefficients af^ = af\u) are 
straightforward to determine using this linear TL algebra realization of the boundary diagrams. 



21 



Simplifying the boundary operator by implementing the condition on the link states yields the 
following proposition. 

Proposition: Implementing the restriction on link states, the boundary operator (A. 5) sim- 
plifies to: 

k (p) (ii t\ _ j + s ( 2u ) p (p) ( a s ^ 

The generalized Temperley-Lieb projectors are 

p-2 k 

P (p) = ^(-l) fc s ((p - k - l)\)e { k \ e ( Q k) = e oei e k = l[e 3 (A.9) 

k=0 j=0 

where sUp + 1)A) = U p ((3/2) are Chebyshev polynomials of the second kind of order p with 
/3 = 2cosA. 

Proof: For one column, that is p = 2, the result follows immediately from the identity 

Xj(u - v)Xj(u + v) = s(u - v - X)s(u + v-X)I + s{2u)e j (A.10) 
with v = £i. For p > 2, we proceed by induction in p using 

^ P+1)(U '° = v^(u% X ° iU " ^ + ^ ( A - H ) 



with 



+ (A.12) 



s(2u) 



Explicitly, 

= s(uH)s(u-QX (u-QX (u+Q + s(2u)X (u-QPi p) X (uH P ) 
= s(u+Os(u-Z P )[s(u+£p-iWu-£ p +i)I + s(2u)e ] (A.14) 
+ s(2u)s(u-Qs(u+^ l )[P^ +1) -U^ 1 (/3/2)e } + s{2u)s{u-Qs{u+QU p _ 2 {(3 /2)e 

= s (u + - e P )[s(M+os(M-e P+ i)/+ S (2 M )p (p+l) ] 

where all terms involving eo cancel out. The term proportional to P[ p ^ is discarded since it is 
killed by the restriction on link states. 

The first few generalized TL projectors are 

Pj 2) = e j} P> 3) = f3e j -e j e j+u Pj 4) = (/3 2 - l)ej- /3eje j+1 + e j e j+ ie j+2 (A. 15) 
These generalized projectors satisfy the properties 

p 0O p( P ) = Up ^{(5/2) pf (A. 16) 

P^eo = 0, e P[ p) = C/ p _ 1 (/3/2)e - P (p+1) , e P[ p) e = U p ^/2)e (A.17) 

22 



The latter properties were used in the proof of the proposition. Although we call the Pj 
projectors, strictly speaking they can only be properly normalized in the case f/ p _i(/3/2) ^ 0, 
that is, p not a multiple of p' . 

Corollary: The expansion coefficients cif^ appearing in (A. 7) are given by 
(p) = 1 ( P ) _ (-iy- 1 s{(p-t)X)s(2u) _ {-\)^U p _ i _ 1 {p/2)s{2u) 

(A.18) 

This follows since opening up the planar diagrams on the right side of (A. 7) for i = k + 1 = 
1, . . . , p — 1 gives the ordered TL products . This also establishes that Kq p '(0, £) = /. 



B Proof of the Inversion Identity 

In this appendix, we prove the inversion identity (4.1) in the planar algebra. Diagrammatically, 
the left side of (4.1) corresponds to 




C 



c 





row 4 






row 3 






row 2 






row 1 





) c 



) c 



— u 




—u 


■v-Zp-i 




-u-£i 


A-Ki 




A-Ki 








X—u 




X—u 


^-£p-2 






u 




u 













(B.2) 



Here and in the following, all face operators have their orientations marked in the bottom-left 
corner. It is also noted that, according to the discussion following the push-through property 
(2.14), we do not need to worry about the restriction rule along horizontal interface between 
D(u) and D(u + A). 

Now, following the proof of the inversion identity in [9], we use the local inversion relation 
to insert a pair of faces with spectral parameters ±2w between rows 2 and 3. We apply the 
push-through property of the YBE repeatedly to bring these faces to the left and right edges. 
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Expanding these faces using (2.3) then transforms the planar diagram (B.2) into 



sin 2u 

cos 2 2u 






row 4 






row 2 






row 3 






row 1 






+ 



sin 2u cos 2u 
cos 2 2u 






row 4 






row 2 






row 3 






row 1 





cos 2u sin 2u 
cos 2 2u 



c 





row 4 






row 2 






row 3 






row 1 






+ 



cos 2 2u 
cos 2 2u 



c 



c 





row 4 






row 2 






row 3 






row 1 





(B.3) 



) 



We now eliminate the last three diagrams by examining the consequences of having a half- 
arc connecting the two left edges (or two right edges) of a two-column appearing in the two 
lower or two upper rows. Expanding in terms of connectivities, we thus use that 



V. 77 



sin rj 



;b.4) 



and 



D= — sin 2 ?] — A — J 



(B.5) 



imply that such a half-arc will propagate. In all three diagrams under consideration (B.3), this 
will eventually lead to a closed loop thereby yielding a vanishing contribution. Also due to the 
propagating property of the half-arcs, the surviving diagram on the right side of (B.3) may be 
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expanded as 



sm 4N u JJjsin(u — sin(u + £,-)] 
i=i 



p-i 



+ 



4JV 
COS U 



\\[cos{u - Q cos(m + Q] 



3=1 




+ 




C7 



zl 



[7 



or 



or 




(B.6) 



In the last diagram, there are four possible connectivities for each four-column configuration 
and they come with the corresponding weights. As demonstrated in the proof of the inversion 
identity in [9], the last diagram in (B.6) is in fact proportional to the identity diagram, as are 
the first two diagrams. Combining the weights following from the various steps in the above 
analysis readily produces (4.3) thus completing the proof. 



C Double-Columns and g-Narayana Numbers 

Definition: A single-column configuration of height M consists of M sites arranged as a 
column. The notion of occupants, signature, weight and corresponding monomial are defined 
as for double-column configurations as in Section 4-3. 

Definition: The polynomial / \ is defined as the sum of the monomials associated to the 

distinct single-column configurations of height M with exactly m occupants. If any of the 
inequalities < m < M is violated, we set ( ) =0. 



Proposition: The polynomial 



M • 
m , 



is given by 



M 



m 



E 



q 



;m(m+l) 



i<h<h<-<im<M 



M 
m 



(c.i; 



Proof: A single-column configuration of height M with exactly m occupants is characterized 
by its signature (j m , . . . ,jx) where 1 < j\ < ji < . . . < j m < M . Summing over all distinct 
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configurations thus amounts to summing over all the associated signatures. This gives rise to the 
multiple sum. The unique column with only the m lowest sites occupied has weight |m(m+ 1), 
while the g-binomial M by construction encapsulates the distribution of excitations of this 

L Jtj 

configuration. This yields the second expression. 

Definition: D M denotes the set of distinct double-column configurations of height M. Simi- 
larly, A M denotes the set of distinct admissible double-column configurations of height M. D^ n 
is the subset of D M consisting of the double-column configurations with exactly m and n occu- 
pants in the left and right columns, respectively. Similarly, A^ n is the subset of A consisting 
of the (admissible) double-column configurations with exactly m and n occupants in the left and 
right columns, respectively. 

Lemma: The polynomial defined as the sum of the monomials associated to the elements of 
Dmn i s given by 



/M\/M\ i m(m+1 ) + i n(n+ i) 
\m /o\ n 1 



M 
m 



M 
n 



(C.2) 



Proof: The distributions of occupancies in the two columns of an element of D^ n are 
independent. The polynomial is thus given by the product of the polynomials associated to 
the individual columns. 



Theorem: For < m < n < M, there is a weight-preserving bijection between D^ n \ A M 
ana u m _ l n+1 . 

Proof: For m = 0, D^ n = A^ n so D^ n \A^ n = = D^_ x n+1 rendering the bijection trivial. 
A non-admissible configuration with m > occupants in the left column violates the inequality 
Lj < Rj for at least one j where 1 < j < m. Let jo denote the minimum such j. We now 
introduce the map 

V^n : (L,R) >-> ((Ri,...,Rj -x,Lj +x,...,L m ),(Li,...,Lj ,Rj ,...,Rn)) (C.3) 

defined for every signature (L,R) corresponding to an element of D^ n \A^ n . Since -R JO -i > 
Lj _i > Lj 0+ i and Lj > Rj , the image of (L, R) is the signature of an element of -D^f-i n+i> so 

,n M ■ n M \ a m — \ n M (c a\ 

This map is clearly injective. To show that it is surjective, we consider the signature (C, 1Z) of 
a given but arbitrary element of -D^f_ x n+l . Since m < n, there exists io, 1 < io < m, such that 
Li > IZi for i < i while £j < 7Zi . We now introduce the map 

<fim-l,n+l '■ (A ^ (C^l> • • • j 'R-ioi A ) • • • j (£l, • • • , A -l> T^ia+U • • • ) T^n+l)) (C5) 

defined for every signature corresponding to an element of D^_ ln+l . As above, it follows that 
the image corresponds to an element of D^ n . To see that it is non-admissible, we note that 
the image in (C.5), here denoted by (L',R'), satisfies L\ < R[ for i < i , while L' io > R' iQ . The 
map 0^-1 n+i * s thus the inverse of the map n and the bijection has been established. The 
bijection (C.3) is weight preserving since the sum of the signature entries of the image of (L, R) 
equals the sum of the signature entries of (L, R) itself. 

Remark: The constraint m < n is important for the existence of the bijection. For m > n, 
A^ n = and -Dm,n\^m,n = ^mn which is generally not isomorphic to D^_ ln+1 . The proof 
above applies for m < n, but breaks down for m > n since the existence of io is no longer 
ensured. 



m.n 
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Remark: There is a very simple configurational realization of this bijection. Viewed from 
above, one identifies the location, jo, of the first violation of admissibility in the non-admissible 
double-column configuration with signature (L,R). One then divides the configuration into two 
with a horizontal cut just below the height where the violation occurred, i.e. just below the 
height Lj (which for Lj = 1 would be just below the entire double-column configuration). One 
finally interchanges the parts of the two columns appearing above the cut while the parts below 
the cut stay put. 

Theorem: The (generalized) q-Narayana numbers defined as the sum of the monomials 
associated to the elements of A^ n is given by 



M 



m, n 



|m(m+l)+|n(n+l) 



~ M~ 




~M~ 




M 




M 




m 


q 


n 


1 


m — 1 


q 


n + 1 


.) 



(C.6) 



M 



can 



Proof: Since D^ n is the disjoint union of A^ n and D^ n \A^ n , the polynomial ^ 

be written as the sum of the monomials associated to the elements of D^ n minus the sum of 
the monomials associated to the elements of D^ n \A^ n . The expression (C.6) now follows from 
the previous Theorem and (C.2). 



D Proofs of ^-Catalan Combinatorial Identities 

The finitized partition functions, for general (r, s) boundary conditions, in (4.34) and (4.35) are 
decompositions into irreducible blocks given by g-Catalan polynomials. In this appendix, we 
prove combinatorial identities to simplify these decompositions to the explicit finitized characters 
given by (4.42). This is achieved by separately expanding out the left and right sides of the 
identities as sums of a single g-binomial. We note that the g-binomials satisfy the elementary 
properties 



1 - q a ~ b + l ) 
(1 - Q b ) 



a 

b-1 



(1 - g°) 

(1 - q a ~ l 



a — 1 
b 



1 - q a ) 

^7) 



a — 1 
b-1 



(d.i; 



and use the notation 



|(iV-2r + s) 



p = 2r 



|(JV-2r + s + l), p = 2r-l 



(D.2) 



We consider four cases depending on the parities of p and s. 
D.I p even, s odd 

We telescope the partition function (4.42) as follows 



3? (9) 



g-4 q^ 2r ~ s)2 



N 
t 



q 



i(2r- s+2) 2 



N 
t-l 



+q 



i(2r+ S -2) 2 



N 
t-s + l 



-q 



+ q 

q 

i(2r+s) 2 



l(2r-s+2) 2 



N 
t-l 



i Wgi 



k=l 



r{2r-s-2+2k) 2 



(I 

N 

t + l-k 



N 
t-s 



±(2r-s+2fc) 2 



t - k 



(D.3) 
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where t = h(N — 2r + s). Using the first part of (D.l) on the first g-binomial, the two terms 
combine 



s i(2r-s+2fc) 2 



ZWfg) = q - l sY— - 

"r,s W « (]_ _ qt+l-k 



q * ^ (1 - 



-Jfc\ 



N 
t - k 

N 
t - k 



(? |(-2r+ S +l-2fc)^ 1 _ q N-t+k^ - (l- q t+l - k ^ 



(q 



l(-2r+s+l-2fe) 



JV+l . 



(D.4) 



Using the definition of the the g-Catalan polynomials in (4.28), the decomposition (4.34) of the 
finitized partition function with p = 2r can be expanded as 



fc=i 



jV-1 2l — a 
2 ' 



^(g) = ^ ? i(2,- S -l+2 fc )(2r- S -3+2*0 
fc=l 



(1 _ g |(2r-s-l+2fc)^ 



iV+l , 



g -I (? §(2r-s+2fc) 2 (? ±(-2r+6-:- 1-2/V. 



l- 2fc )(l-^ (2r - S - 1+2A) ) 



iV + 1 
t + 1 - k 



(D.5) 



k=l 



(1-9 2 ) 



iV + 1 
t+ 1 - k 



Applying the third part of (D.l), we can simplify this expression to obtain the desired result in 
agreement with (D.4) 



_i • q \{2r-s+2kf 



k=l ^ q ' 
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N 
t - k 



( ? |(-2r+ S +l-2fc) _ 1 ^ 1 _ q 
(1-9 2 ) 



(D.6) 



D.2 p odd, s odd 

We telescope the partition function (4.42) as follows 
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(D.7) 



where t = |(JV — 2r + s + 1). Using the first part of (D.l) on the first g-binomial, the two terms 
combine 



|(2r-s+2fc) 2 r 
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it 
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(D.8) 
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Using the definition of the the g-Catalan polynomials in (4.28), the decomposition (4.34) of the 
finitized partition function with p = 2r — 1 can be expanded as 
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(D.9) 



Applying the first part of (D.l) to the second g-binomial, we can combine the two terms to 
obtain the desired result in agreement with (D.8) 
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(D.10) 



D.3 p even, s even 

We telescope the partition function (4.42) as follows 
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(D.ll) 



where t = -(AT — 2r + s). Using the first part of (D.l) twice on the first g-binomial, the two 



terms combine 
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Using the definition of the the g-Catalan polynomials in (4.28), the decomposition (4.35) of the 
finitized partition function with p = 2r can be expanded as 
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(D.13) 



Applying the third and first parts of (D.l), we can simplify this expression to obtain the desired 
result in agreement with (D.12) 
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D.4 p odd, s even 

We telescope the partition function (4.42) as follows 
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where t = |(AT — 2r + s + 1). Using the first part of (D.l) twice on the first g-binomial, the two 
terms combine 
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(D.16) 
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Using the definition of the the g-Catalan polynomials in (4.28), the decomposition (4.35) of the 
finitized partition function with p = 2r — 1 can be expanded as 
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Applying the first part of (D.l) to the first g-binomial, combining the two terms and then 
applying the identity again, we obtain the desired result in agreement with (D.16) 
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